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Heavy quark medium induced radiative energy loss is derived to all orders in opacity, (L/Xg)". The 
0^ . analytic expression generalizes the GLV opacity expansion for massless quanta to heavy quarks with mass 

M in a QCD plasma with a gluon dispersion characterized by an asymptotic plasmon mass, rUg = gT/\/2. 
Remarkably, we find that the general result is obtained by simply shifting all frequencies in the GLV series 
by {rn^ + x'^ M'^) / {2xE). Numerical evaluation of the first order in opacity energy loss shows that both charm 
and bottom energy losses are much closer to the incoherent radiation limit than light partons in nuclear 
collisions at both RHIC and LHC energies. However, the radiation lengths of heavy quarks remain large 
■ compared to nuclear dimensions and hence high pr heavy quark production is volume rather than surface 

' dominated. 

o 



1 Introduction 

The discovery ^-[51 of a factor of 4 ~ 5 suppression of high p_L ~ 5 — 10 GeV hadrons produced in central 
O , Au + Au at the Relativistic Heavy Ion CoUider (RHIC) has been interpreted as evidence for jet quenching of 

light quark and gluon jets 0-^1 ■ Jet quenching was predicted [201 to occur due to radiative energy loss of high 
energy partons that propagate through ultra-dense QCD matter. The observed quenching pattern is interesting 
not only as a test of QCD multiple scattering theory but also because it provides a novel tomographic tool 
that can be used to map the evolution of the quark gluon plasma plasma (QGP) produced in ultra-relativistic 
• nuclear collisions. Medium induced radiation arises from higher twist final state interactions and depends on 

5^ \ the integrated optical thickness or gluonic opacity L/\g of the QCD medium. At high pr it dominates over the 

energy loss due to elastic scattering |21| . 

Recent data from RHIC D + Au reactions on the other hand, show the absence of jet quenching 

in such light ion beam reactions. These data are in accord with predictions based on jet tomography in Refs. 
|10[ nil [T^ because in light ion interactions no extended dense QCD medium is produced in the final state. 
This important D + Au control experiment therefore strengthens the interpretation of the observed high pT 
quenching pattern in central Au + Au as due to final state energy loss of jets produced in dense QGP matter. 

The first estimates for heavy quark energy loss |ini-|2HI proposed that similar quenching may occur for 
charm jets as for light partons. However, in ref. [221 it was pointed out that the heavy quark mass leads to 
a kinematical "dead cone" effect for 9 < M/E that reduces significantly the induced radiative energy loss of 
heavy quarks. Numerical estimates indicated that the quenching of charm quarks may be about a half that of 
light quarks. Experimentally, the first PHENIX data PO] on "prompt" single electron production in Au + Au 
collisions at ^/s — 130 AGeV provided a first rough look at heavy quark transverse momentum distributions at 
RHIC. Remarkably, no indication for a QCD medium effect was found within the admittedly large experimental 
errors (see also ref. (31)'). However, in the near future, data with much higher statistics and wider pT range will 
become accessible. 
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We concentrate in this paper on the theory of heavy quark energy loss, extending our previous works I |32| 
and II ''6'6K Their quenching pattern and correlations should provide independent complementary tests of QGP 
production. In addition, heavy quark observables can be used to test different approximations in the emerging 
theories of heavy quark radiative energy loss QCD [201-1221 i The main new feature that heavy quark 

probes introduce is a controlled reduction of the radiated gluon formation times due to finite mass kinematic 
effects. As the mass increases, the increase of phase shifts leads to a reduction [221 destructive Landau- 
Pomeanchuk-Migdal (LPM) interference effects that were found to be so important for light quark or gluon jet 
energy loss "^-'IQ'. In addition, the kinematic reduction of propagator amplitudes leads to an overall decrease 
of the magnitude of heavy energy loss relative to light jet energy loss even in the incoherent Gunion-Bertsch 35j 
(GB) limit, as we demonstrate quantitatively below. Therefore, the predictable high px open charm and open 
bottom observables at RHIC and LHC will provide interesting new control tests of QCD dynamics in nuclear 
collisions. 

The main goal of this paper is (1) to generalize the GLV opacity series to include massive quark kinematic 
effects and (2) to take into account the Ter-Mikayelian plasmon effects for gluons as described in Ref. [33|. The 
competition between these two medium effects was first discussed in our paper I |32j . We showed that the 
apparent null effect observed for heavy quarks via single electrons could in part be due to a reduction of the 
leading 0{x'^) order radiation associated with the initial hard production process. In this sense, there are two 
opposing medium effects: (1) at 0{x'^), the Ter-Mikayelian plasmon dispersion effect that reduces the associated 
hard radiation energy loss [33| and (2) at 0(x"-^), the induced radiative energy loss that increases the total 
radiative energy loss albeit with a reduced efficiency due to the dead cone effect ^S*. The detailed derivation of 
the QCD Ter-Mikayelian effect was presented in paper II |33| . An important conclusion from II, which we rely 
on in the derivation below, was that the effects due to the plasmon dispersion relation can be well approximated 
for high px jets ignoring the longitudinal modes and applying the asymptotic (short wavelength transverse) 
plasmon mass. That asymptotic mass, nig = /i/\/2 is somewhat larger than the long wavelength plasmon mass 
n/y/S, where /x « gT GeV is chromoelectric Debye screening mass. 

We concentrate here as in ITS' on the case of high pT eikonal jets produced inside a finite evolving QGP at 
some initial point (to-, zcXq). This is in contrast to the Gunion-Bertsch problem 35.^ with to = — oo, where the 
jet is prepared as a beam in the remote past. We model the interactions in the QGP as in I17j via random 
color screened potentials 



where x„ is the location of the n**^ scattering center and w(qn) = 47ras/(q^ + M^)- The elastic cross section 
between the jet and target partons in the GW model is 

d<JeiiR,T) ^ CrC2{T) \v{q)\^ 
d2q dA (27r)2 ■ 

For consistency with GLV we use the same notation throughout. Transverse 2D vectors are denoted 
as bold p, 3D vectors as vectors p = {pz,p), and four vectors hy p — {p°,p) — + p^ ,p^ — p^ ,p]. The color 
exchange bookkeeping with the target parton n is handled by an appropriate SU{Nc) generator IKj, Ta{n), 
in the dn dimensional representation of the target (TrTa(n) — and Tr (Ta(«)Tb(j)) = 5ij5abC2{i)di/ dA)- We 
assume that all target partons are in the same dx dimensional representation with Casimir C2(T).) We denote 
the generators in the du dimensional representation corresponding to the jet hy a = ta with aa = CrI. For 
heavy quarks in SU{3), Cr — 4/3 while Ca = 3. The elastic cross section of target parton i with the jet is 
therefore proportional to the product of Casimirs, CrC2{T). As in the analytic results derived below with 
the reaction operator approach do not depend on the actual form of v, but the Yukawa form will be used for 
convenience in numerical estimates. 

The sections are organized as follows: In section 2, we review the zeroth order in opacity but in medium 
associated radiation for a massive quark jets. At zeroth order only the hard initial vertex acts as a source 
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for gluon radiation, but in a medium the dynamical gluon mass, \J^'^i{k) — fc^ ~ jji^ jg taken into account 

|33| . In section 3, the first order in opacity heavy quark energy loss is computed. We improve Eq. (2) in 
paper I [32] ( see Eq. (|12|l below) by incorporating more carefully the kinematic constraints so that the GLV 
massless M ^ limit can be recovered. The detailed diagrammatic computations are recorded in Appendices 
A-G closely following the derivation in GLV [TS] but pointing out the different phases and residues that arise 
when M and rUg are not zero. Section 4, presents the generalization of GLV to all orders in opacity for finite 
masses. Our main result fEa. (|17|l ). which follows from the reaction operator method using the results from the 
appendices, is that the analytic structure of the small-x induced radiation distribution in the heavy quark case 
is identical to that of the massless case but with a single universal energy shift Aw = (m^ + M"^) / {2xE) . The 
energy shift is independent of the momentum transfers qi and the plasmon transverse momentum k. In section 
5 we give numerical estimates for the the first order energy loss of charm and bottom quarks comparing RHIC 
and LHC conditions. We show quantitatively how close our (partially coherent) numerical results are to the 
incoherent limit. Conclusions are presented in section 6. 



2 The one gluon associated radiation 

In order to compute the associated radiative energy loss when the hard process is embedded in a dielectric 
medium, we need to compute the squared amplitude of Feynman diagram, Af°^^, which represents the source 
J that produces an off-shell jet with momentum p' and subsequently radiates a gluon obeying the dispersion 
relation, uj{k), of the medium with momentum k. The jet emerges with momentum p and mass M. Since our 
focus is on heavy quarks, we neglect the thermal shifts of the heavy quark mass here. As noted before, the 
detailed computation of the associated energy loss was done in ^33', where it was shown that gluons in the 
medium can be approximated as massive transverse plasmons with mass rUg w TOoo = fJ-fV^. 

If we assume that pz and kz are large enough, such that M'^/pl ^ 1 and m?g/k1 <^ 1 are satisfied, then we 
can write p, fc, and the transverse polarization e in terms of light cone components: 



k^[xE+,k-M ,e(fc) = [0,2^,e] , p = [(1 - x)E+, p", p] . (3) 

Soft radiation is defined as a; ^ 1 so that, for example, p+ ^ fc+ and p^ = (p^ + A/^)/(l — x)E^ <^ k^ = 
(k^ + m?g)/xE^ . Here, M is the mass of the quark, and mg is the mass of the glue of energy uo ~ xE^ /2. We 
also adopt the same shorthand notation as in for energy differences: 



k2 (k-q.)^ (k-q.-q,-)^ 



In the soft eikonal kinematics that we consider 

r.+ ,+ ™9 (p + k)2+M2 

^+ » fc+ » w(,...j) + ^ » ^ . (5) 

The hard jet radiation amplitude to emit a transverse plasmon with momentum, polarization, and color 
(fc, e, c) without final state interactions is 



M,",rf = zJ(p + fc)e*(''+'=)"«(z5,)(2p + fc)^e^(fc)zAMb + A)c 

ek 



J(p)e^''-»(-2zg,) 



A:2 + m2 + M2a;2 



e''^«^° c , 



(6) 
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where a;o(0, zq, 0) is the jet production point inside the plasma. We assume, as in "16 , that J varies slowly 
with p, and neglect high x spin effects. Eq.® represents the mass corrections to the Eg. (35) in 16 . In soft 
radiation approximation the spectrum can be extracted as 



2£;(27r)3 2w(27r^3 
where (with dfl = 3 dimensional representation quarks) 



Eqs. {THHJ lead to the finite mass generalization of the small x invariant DGLAP radiation spectrum 

(9) 



MU \fLso^fLs - d'Njd'Nl,^^ , (7) 



dN^°^ dA^^°^ Cnas 

LJ ^ K. X- ~ 



d^k dxd^k^ (k2+m2+a;2M2)2 
Eq. (jHJ clearly shows the depletion of radiation in the "dead cone" at angles 

e kOc^ ^m?g + x^M^/{xE) 
generalized to take into account also the Ter-Mikayelian dielectric dispersion in a QGP. 

3 First order radiative energy loss 

The first order in opacity energy loss can be computed from formula: 



d^k 



dr ' "''(27r)32pO (27r)32tj 



(10) 



where d^Nj is given by Eq. © and dx is the dimension of the target color representation {— 8 for a pure 
gluon plasma). Afi is sum of all diagrams with one scattering center and M2 is sum of all diagrams with two 
scattering centers in the contact limit. We compute those diagrams using the same assumptions as in [16] . as 
reviewed in Appendix A. The detailed evaluation of the amplitudes is presented in Appendices B-F. The results 
are combined in appendix G, to give the small x differential energy loss 



d^ ^ CnasL^f dk^ r d\^ 

dx TT A 7 k2 + m2 + M2a;2 J n {qj + ^2)2 

k ■ qi(k - qi)2 + (m2 + M2a;2) qi ■ (qi - k) 
X (4|^)2 + ((k_q^)2 + M2a;2 + m2)2 ' ^ ' 

where |qi| is the magnitude of the transverse momentum transfer between a target parton and a jet, |k| is the 
magnitude of the transverse momentum of the radiated gluon, and Xg is the mean free path of the gluon. The 
simple analytic form of the destructive interference factor involving {AEx/L)'^ arises for an assumed exponential 
distribution exp(— Az/L)/L of the distance between the jet production and target rescattering center. Note 
that in the massless limit, M = nig — 0, Eq. (|ll|l reduces to Eq. (125) in jl6|. 
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Numerical evaluation of Ea llll shows that to a good approximation we can ignore the finite kinematic bounds 
on |qi| < V6-ET for E' > 10 GeV charm jets. With this simplification, we can substitute q ^ qi — k and 
integrate over the azimuthal angle to arrive at 



AE^^l = d.Er 



Jo Jo + (q' + M^x^ + mlY 



((|k|-|q|)2 + ;,»)-VJ((|k| + |q|)» + ,. 2)3/2 _ 

If we further neglect the finite kinematic boundaries on the gluon transverse momentum, |k| < xE, then in 
the massless M = nig — limit, we recover the approximate asymptotic Eq. (127) of The integral in 
Eq. H12I) can be performed analytically if as does not run, but the cumbersome result is not instructive. 

We note that Eq. lO differs from the Eg. (2) in 33 . The A£'|^^ in [S] was obtained by assuming that 
energy of the glue is much larger than momentum transfer qi, i.e xE ± |qi| « xE. Eq. H12|l avoids this 
approximation by incorporating exact kinematics at the price of a more complicated expression. Numerically, 
this improvement however, does not change significantly the estimated charm energy loss presented in |33) as 
shown in section 5. However, theoretically, Eq. (|12|l is prefered since the Af — s- limit is only correctly recovered 
when exact kinematics are enforced. 

The main qualitative effect of increasing M is to reduce the relevance of the inverse formation time factor 
[AEx/LY in the denominator of the integrand. Formally, by setting this factor to zero recovers the incoherent 
limit of induced radiation. This corresponds to the QCD analog of the QED Bcthe-Hcitler limit with AEqb 
agEL/Xg modulo a logarithmic factor. In QCD this includes the isolated scattering Gunion-Bertsch (GB) 
radiation as well as elastic scattering of the associated radiation from zcroth order in opacity as we discuss in 
the following section. 

The incoherent, short formation time limit, generalized here to include the plasmon asymptotic mass, is 
given by 

AEin Cpas L f\ f 2qV^dq' 

— ax 



E XgJo J (q2 + M2a;2 + m2)2 

dk"^ e{xE - \k\) j 2 .,2 ^ k"^ - M^x^ - ml'] 

((|k| - |q|)2 + Ai2)3/2((|k| + |q|)2 + ^2)3/2 ^ ^ ^ ^ ' k^ + APx^ + j ' 

Since we found in Ref. and show in more detail in section 5 that the heavy quark energy loss is close 
to the incoherent limit, it is also useful to define the effective gluon radiation length, Lrad, via 

dAE E , , 

= (14) 

dL LradiM,mg,E,L) 

This length is of relevance to answer the question whether the jets observed in a given kinematic window 
are dominated by surface or volume emission. 

4 Higher orders in opacity energy loss 

4.1 Heavy quark generalization of GLV to all orders 

From the results of the calculations reported in the appendices, we find that the finite masses modify the 
effective radiation amplitudes and phase factors in a remarkably simple and universal way. The phase factors in 
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the massless case studied in GLV modulate the amplitudes by factors such as Exp{i a;(„...„-) Az™). The 
energy differences W(m - n) from Eq.l^J are inverse formation times and Az^ are distances between scattering 
centers. With finite M and nig these energy differences are found from Eqs. (|^ HSl ICT EOl EHl EZI [Z2Il to be 
simply shifted by a independent term, 



_ (k-q^ q„)^ _ m|+M^ 

- 1^ * "(m.-,«) = ^(m,-,n) + U^j 

In addition, the kinematic current amplitudes appearing in those equations are simply modified versions of 
the Hard, Gluon-Bertsch and Cascade terms in GLV, which for finite masses are now 



k - ^ (k-Qn -Qi, qz„J 

k2 + ml + A'Px^ ' (n^2■■■.™) (k - q,^ - q,^ q,^)2 + m2 + A'Px^ ' 

In summary, the computations in Appendices B-F show that the diagrams for the finite masses case can 
be obtained from corresponding massless equations in jl5) . by simply replacing the terms from Eqs. (58, 106) 
in by the modified terms defined in Eqs. ifTBlIT^ . 

Therefore, the recursive GLV reaction operator formalism carries over to the massive case with the simple 
replacements above. The complete arbitrary order in opacity induced radiation distribution can be obtained by 
generalizing Eq.(113) of GLV as follows: 



Bj = H — Ci , 



dx (Pk 



-2C(i,...,„) • B 



cos ^ 17(fc^...^„)Azfc -cos n(k,---,n)^Zk 



, k=2 / \ k=l 



(17) 

where ^2 = is understood and |ui(qi)p is defined as the normalized distribution of momentum transfers 
from i^^ scattering center ^Hl- For the Yukawa screened interactions (Eq. Q), the differential gluon cross 
section in the local density approximation is 

a(..q.) ^ a.(z.)|..(qOI = = . {^^ + .{z.W ' ^''^ 

where fJ.{zi) is the local Debye mass that may vary if the system expands. 

Note that, as in ^B] the Eq. (|17|) is not restricted to uncorrelated geometries. Also it allows the inclusion 
of finite kinematic boundaries on the q^ as well as different functional forms of the gluon elastic cross sections 
along the eikonal path. The n}^ order in opacity energy loss spectrum can be obtained from Eq. 1)17(1 via 

= / 0? k— — xE (19) 



dx J dx d'^'k 

From the Eqs. 1(171 119|l we can obtain the first order in opacity energy loss by setting n = 1. The result 
obtained after averaging over exp(— Azi/i)/L is the same as Eq. ((84(1 in Appendix G, which leads to Eq. ((12(1 
computed in previous section. 
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In order to make the averaging over the target coordinates more explicit, assume an uncorrelated geometry 
and let p{z) = p{z,T = z) denote the target density along the path of the jet. For Bjorken expansion p(z) = 
6{R — z)poTo/ z for example. In the local density approximation, the screening mass, may also depend on 

the proper time (r = z). Then (J(.i{z) and q)p may vary along the jet path as well. The average over over 
the target centers can be made explicit by replacing the opacity factor in (|17|l by 



dzip{zi) ■■■ dZnp{Zn) X 

p n 

J n ('^'q^ H^^,^^) - <Jel{z)S^{(^^)]) (20) 



4.2 The Incoherent, Short Formation Time Limit 

The incoherent limit of Ea. (|17|l is obtained formally by taking the large r2(fc^....„) Azfe ^ 1 limit in which all 
cos(-) factors average to zero except for one term for m = 1 in which the cos(O) = 1. Hence the incoherent limit 
of the n*'' order in opacity induced radiation is 



dx (i^k 

For the n — 1 contribution 
r(i) 



^""^ ^-^j I d'<iiv'M(Bi„m + cU-W] , (22) 



where B(o)(i) = H — C(i) is the finite mass generalization of the incoherent GB radiation amplitude including 
the asymptotic plasmon dispersion. Note that (see Wiedemann 23|) the — term corresponds to the first order 
unitarity correction of the full Glauber associated zeroth order in opacity contribution e^^^^^'dN'^^^ /dxd'^h. In 
addition elastic rescattcring of the radiated gluon leads to the cascade contribution C^j^j which further broadens 
the (already very broad) transverse momentum distribution of the associated radiation. The incoherent induced 
radiation from the rescattcring of the quark jet in the medium is the Gunion-Bertsch term B^q^^j^-j. There are 
L/Xg such induced contributions in the completely incoherent limit. 

In general we expect that with the inclusion of masses in Eq.lUnj, the increase of the r2(m ..._„) Az„i arguments 
of the interference cosines in Ea. H17|) will drive the radiation distribution closer to the incoherent limit. In the 
next section we will evaluate numerically how close is charm and bottom quark induced radiation to their 
respective incoherent limits. 

Similar to the massless case studied in GLV I^, we expect that in the case of finite mass, higher orders 
in opacity contributions to the net induced energy loss will converge rapidly for moderate L/Xg of practical 
interest. The dominance of the first order contribution to the transverse momentum integrated energy loss is 
more readily S66n in tliG incoliGrent limit Rbove. CliRiige variciblGS k — >■ k' = k — Q2 — ... — q,^ 

and integrate 

over k'. In the lj ^ p limit, we can approximately ignore the change in the |k'| < lo kinematic limit. In that 
case 



•^^k C(i ..._„)• B(2,...,n)(i,...,„) « / d^k' C(i) • B(o)(i) . (23) 
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However, the integrated energy loss at order n > 2 has additional integrations 

P n 



i=2 

that tend to kill the whole contribution modulo tiny edge effects due to Ea. l23|) not being exact. Therefore, 
especially as we approach the incoherent limit the main order contribution to the medium induced energy loss 
is dominated by the first order in opacity. 

In order to see how this works in more detail consider the case of uncorrelated geometries where we can use 
Ea. (|2()|l and Fourier techniques to sum the whole opacity series in this incoherent limit. Note first that the term 
proportional to (5^(qi) vanishes due to B(2,- -,n)(i,- -,n) = when qi = 0. It is convenient therefore to define the 
accumulated Q — 1* momentum transfer distribution by inserting a factor 

1 = I d'QS'iQ - q,) ^ I Jd'h e»Q-^ n e-'<^--b 

Define an effective dipole cross section 

and a corresponding dipole opacity 

Xd{z,h) = p{z.,)<Td{z,h). 

In terms of these quantities, 



dx d^k 



(2^)2 



dzid^qip(zi)cr(zi,qi) 
d^be^'^-b / dz2{~l)xd{z2,b) •••/ dzr.{~l)xd{zn,b) X 



x2Ck-Q-qj • (^Ck-Q-qj — Ck-Q 

dzid^qip[zi)a{zi,qi) 
^'Q /'^2K„.Q b (-1) 



(2^)2 

x2Ck-Q-qj • (^Ck-Q-q^ - Ck-Q^ (24) 

We can now sum all orders from n = 1 , oo into a closed form 

dNin CrUs 



dx d^k 



/ (izid^qip(zi)cr(zi,qi) 



(2^)2 



■2i iQ-b ~ dz'p{z')(Tdiz',b) 



d^h 



e 



x2Ck-Q-q^ ■ (^Ck-Q-qj - Ck-qj (25) 

We can now clearly see how (|25|l reduces to the first order result when integrated over k. If we change 
variables to k' = k — Q and assume that k' < uj remains approximately valid, then the Q integration reduces 
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to a (5^(b) which in turn converts ad{z,b) — > ad{z,0) = 0. The whole second line reduces then to a unit factor 
so that 



dNi„ 
dx 



-10 



dzd^qp{z)aiz,q) / d^k 2 Ck_q^ • Ck-q, -H 



(26) 



However, in the same approximation that led to H26(l from H25(l . we can again shift variables to k' = k — qi 
and again neglect the change in the \k'\ < ut bound to cancel the & — contribution. This collapses the full 
result to the simple incoherent Gunion-Bertsch limit 



dN,„ 
dx 



dNcB CrUs 



dx 



dzd'^qp{z)a{z,q) / d^k B^o^j^) 



(27) 



In order to assess the accuracy of the approximations involved in the approximate handling of the kinematic k 
bounds, we have to evaluate numerically the difference between dNm and dNsG- 

5 Numerical estimates 



5.1 Heavy quark energy loss at RHIC 

The numerical results for the first order induced radiative energy loss are shown on Fig.l for charm and bottom 
quarks. From the analysis of light quark quenching in Au + Au at 130 GeV jHOl effective static plasma opacity 
L/A is in the range 3 — 4. On Fig.l we fix opacity to L/\ = 4, and we look at the 1** order in opacity fractional 
energy loss as a function of initial energy of the quark. We assume here that ag — 0.3, p = 0.5 GeV, and A = 1 
fm for the plasma parameters. Since finite parton masses shield coUinear k — s- singularity |33) . our numerical 
computations are performed with zero momentum cutoff. We see that for heavy quarks, in the energy range 
E ^ 5—15 GeV, the Ter-Mikayelian effect reduces the induced energy loss in this extension of the GLV approach 
somewhat more than in the BDMS approximation |29j . However, on an absolute scale, this only corresponds 
to a change of S{AE'^^^ /E) < 0.05, which is negligible. Note that with both dead cone and plasmon mass 
reduction, there remains a sizeable induced energy loss fraction AE^^^ / E w 0.15 for charm quarks while only 
about half that is predicted for b ottom. 




15 20 
E [GeV] 

FIG 1. The 1"' order in opacity fractional energy loss for heavy quarks with (solid curves) and without 
(dashed curves) Ter-Mikayelian effect approximated by a constant nig — /x/\/2 for transverse modes only. 
Upper curves correspond to charm, and lower to bottom quarks as a function of their energy in a plasma 
characterized by as = 0.3, n — 0.5 GeV, and L = 4A = 4 fm. 
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However, we want to emphasize that, though Ter-Mikayehan effect is not important for the medium induced 
energy loss of the massive quarks, it is very important for massless quark case. Note that excluding Ter- 
Mikayelian effect means setting both plasmon mass and momentum cutoff to zero. Since in the massless parton 
case, k ^ singularities are not naturally regulated, excluding the Ter-Mikayelian effect would lead to infinite 
result of medium induced energy loss. In '16' these divergences are prevented by introducing the finite |k| > /i 
cutoff. It is easy to see from numerical computations that medium induced energy loss for light quarks with 
finite plasmon mass mg = /i/V2 and zero momentum cutoff, is similar to the energy loss obtained with zero 
plasmon mass and finite momentum cutoff |k| > /i. Setting the finite momentum cutoff in the case of 
light partons thus numerically produces similar results as introducing the finite plasmon mass, i.e. including 
the Ter-Mikayelian effect. Therefore, though the Ter-Mikayelian effect was not considered in ^H], the same 
effect was obtained by including the finite momentum cutoff. However, we would like to note that, though 
both approaches give similar numerical results, the advantage of the Ter-Mikayelian effect is that it provides 
the natural regulation of k — > divergences. 

Fig. 2 shows the first order induced radiative energy loss for charm and bottom quarks as a function of 
opacity. We see that the induced contribution increases with L. Note that for heavy (charm and bottom) 
quarks the thickness dependence is closer to the linear Bethe-Heitler like form, L^, than the asymptotic energy 
quadratic form |14l llfij . 




1 2 3 4 5 
L [fm] 

FIG 2. The 1" order in opacity fractional energy loss for a 10 GeV charm and bottom quark is plotted 
versus the effective static thickness L of a plasma characterized hy fi — 0.5 GeV and A = 1 fm. Upper 
curve correspond to charm, and lower to bottom quark. 



In the case of finite kinematic bounds on |qi| < \/6ET, the angular integration in Eq. (|ll|l can still be 
performed, but the result is cumbersome. However, including finite kinematic bound on |qi| is found to be 
unimportant, because it only reduces the induced energy loss fraction AE'-^'^/E in Fig. 1 by less than 10%. 
Therefore, we can safely neglect the finite kinematic boundaries on the momentum transfers. 
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AE 



(1) 
ind 




3 

L[fm] 

FIG 3. Comparison of the incoherent limit (dashed curves) to the 1"* order (solid curves) in opacity 
fractional energy loss for i? = 10 GeV charm quarks via Eq. 11121 under RHIC conditions {as — 0.3, fj, = 
0.5 GeV). The higher dashed and solid curves correspond to nig — 0, while the lower to nig = /^/\/2. 
Arrows point between corresponding incoherent limit and full 1"* order results. 



In Fig. 3 the first order energy loss for charm quarks is compared to the incoherent limit obtained from Ea. H12|l 
by setting AxE/L = 0. This analogue of the QED Bethe-Heitler limit is the one where no destructive LPM 
interference effects occur in QCD. The characteristic of the incoherent limit is the complete linear dependence 
of AE on L. The slope depends on the assumed fj, as well as the assumed dynamical gluon mass. Two 
cases, nig — 0, fi/V^ are shown by the two dashed lines. Comparing to the soHd Hues that include destructive 
interference effects due to finite formation times, we see that the realistic case with finite nig ~ 0.35 GeV is 
remarkably close to the incoherent limit for E — 10 GeV quarks. 

We have also checked numerically the difference between the full first order incoherent limit and the isolated 
Gunion-Bertsch contribution. These should be approximately the same if the elastic rescattering and the hard 
unitarity correction terms cancel as in Ea. (|27|) . We find that the GB energy loss for E=10 GeV charm quarks 
is AEgb/E « 0.35 is significantly larger than the incoherent energy loss AEin/E w 0.22. However by -E = 30 
GeV, the GB limit exceeds the incoherent limit by only 15%. 

Finally, by using the results for associated energy loss from |33j . we can now compare the net (AE^^^^ + 

^^med) energy loss in the medium, with the one in the vacuum defined by fivac = GeV, as shown on Fig. 4. 
As in paper I, we see that even in the absence of a medium (L = 0), a charm quark with energy ~ 10 GeV suffers 
an average energy loss, AeI'^Jc/E « 1/3, due to the sudden change of the color current when it is formed in the 
vacuum. The dielectric plasmon effect reduces this to about AE^^^^/ E « 1/4. This contribution is independent 
of the thickness of the plasma as long as L is not too small. For very small L < 1/mg, the plasmon dispersion 
is smeared out due to the uncertainty principle, and AE^^^^^ must approach AE^ac from below. 
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FIG 4. The fractional energy loss for a 10 GeV charm quark is plotted versus the effective static thickness 
L of a plasma characterized by /i = 0.5 GeV and A = 1 fm. The dashed middle horizontal line corresponds 
to the energy loss in the vacuum taking into account the kinematic dead cone of radiation for heavy 
quarks |29j . The lower horizontal solid line shows our estimate of the reduction of the zeroth order energy 
loss due to the QCD analog of the Ter-Mikayelian effect. The solid curve corresponds to the net energy 
loss, A£i°L + Ai5£' . 



We see that the net energy loss, AiJ^^^ + Ai?|^^, is found to be even smaller than in PSI due to the 
corrections made in Eq. (|12|l . Also, notice that for the effective static opacity medium [3^ with L in the range 
of 3 — 4 fm, /i « 0.5 GeV and A « Ifm, the difference between the medium energy loss (AE'^^^ + AE^^^) and 

naive vacuum value {AEv^c) is small, i.e. between 5% and 10% of the initial energy of the quark. These results, 
therefore, suggest that in addition to the heavy quark dead cone effect, the apparent null effect observed for 
heavy quark energy loss via single electrons may in part be due to a further reduction of the leading order energy 
loss. The different dependence of these effects on the plasma thickness L and on the transport properties, /x 
and A, should make it possible to test this explanation by varying the beam energy. A, and centrality at RHIC 
and LHC energies. 

5.2 Heavy quark energy loss at LHC 

In this section we want to estimate the heavy quark energy loss at LHC. According to some estimates 0, the 
density of plasma partons at LHC should be 3-4 times larger than the one at RHIC. Therefore the results bellow 
are calculated using A and fi based on the estimate given above. 

The numerical results for the first order induced radiative energy loss are shown on Fig. 5 for charm and 
bottom quarks. We assume here that as = 0.3, fj, = 0.7 GeV, A = 0.7 fm and L = 4 fm for the plasma 
parameters. As before, the Ter-Mikayelian effect does not have significant effect on the induced energy loss. On 
the other hand, we see that the induced energy loss fraction is two times larger than before, i.e. AE'^^^ / E w 0.3 
for charm quarks and about half for bottom. 
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5 10 15 20 25 30 
E [GeV] 

FIG 5. The I''* order in opacity fractional energy loss for heavy quarks with (solid curves) and without 
(dashed curves) Ter-Mikayelian effect approximated by a constant nig — ii/\/2 for transverse modes only. 
Upper curves correspond to charm, and lower to bottom quarks as a function of their energy in a plasma 
characterized by as = 0.3, yu = 0.7 GeV, A = 0.7 fm and L = 4 fm. 

Fig. 6 shows the first order induced radiative energy loss for charm and bottom quarks as a function of 
opacity. We see that, again, for bottom quarks the thickness dependence is still closer to the linear Bethe- 
Heitler like form, L^, than the asymptotic energy quadratic form jl4l llfij . However, for charm quark, the 
thickness dependence is somewhere between linear and quadratic form. 




[, I 

1 2 3 4 5 
L [fm] 

FIG 6. The 1"* order in opacity fractional energy loss for a 30 GeV charm and bottom quark is plotted 
versus the effective static thickness L of a plasma characterized by /j, = 0.7 GeV and A = 0.7 fm. Upper 
curve correspond to charm, and lower to bottom quark. 
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FIG 7. Comparison of the incoherent limit (dashed curves) to the 1"* order (solid curves) in opacity 
fractional energy loss for i? = 30 GeV charm quarks via Eg. 1121 under estimated LHC conditions (as = 
0.3, \i = 0.7 GeV). As in Fig. 3, the higher dashed and solid curves correspond to nig — 0, while the lower 



to rug = fi/V2. Arrows point between corresponding incoherent limit and full 1° 



order results. 



Comparing Fig. 7 to Fig. 3, we see that for 30 GeV quarks, the first order result is now further away from the 
incoherent limit. The destructive LPM finite formation effects included in the full 1** order result substantially 
reduce the incoherent energy loss. 

For LHC conditions, we can again, by using the results for associated energy loss from 33 , compare the net 
{AE^^j^ + AE'^g^) energy loss in the medium, with the one in the vacuum defined by ^yac = GeV, as shown 

on Fig. 8. We see that the dielectric plasmon effect reduces the AE^^^^jE form k. 1/3 to w 1/5. As before, this 
contribution is independent of the thickness of the plasma as long as L is not too small. 
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FIG 8. The fractional energy loss for a 30 GeV charm quark is plotted versus the effective static thickness L 
of a plasma characterized hy ^ = 0.7 GeV and A = 0.7 fm. The dashed middle horizontal line corresponds 
to the energy loss in the vacuum taking into account the kinematic dead cone of radiation for heavy 
quarks I29| . The lower horizontal solid line shows our estimate of the reduction of the zeroth order energy 
loss due to the QCD analog of the Ter-Mikayelian effect. The solid curve corresponds to the net energy 
loss, AE^2, + ^Ef,^,- 



We see that the difference between net energy loss, Ae'^I^ + Ai?|^^ and naive vacuum value AE^^c at 
LHC is two times lager than the same difference at RHIC. Therefore, these results suggest a somewhat larger 
suppression of charm quark at LHC compared to RHIC. 



6 Summary 

In this paper wc generalized the GLV radiative energy loss formalism to heavy quarks including also the plasmon 
effects. We have shown that the reaction operator method introduced in jE] can also be applied to the case of 
massive quarks and gluons. Remarkably, a simple mass dependent energy shift, via Eas. l|15ll6|l . was found to 
modify all direct and virtual diagrams. This result, proven in the appendices A-G made it possible to generalize 
the GLV zero mass results to the case of heavy quarks to all orders in opacity. We also derived the incoherent 
limit of the induced radiation spectrum from heavy quarks and found a closed form expression summed to all 
orders in opacity. 

The dependence of the charm first order in opacity energy loss on the jet energy and plasma parameters 
were studied numerically. We showed that the Ter-Mikayelian effect on the induced energy loss is comparable 
but somewhat larger than in (29) . In addition, it was shown shown that the induced contribution increases 
approximately linearly with L, as first reported in |33| . For charm quarks at RHIC the heavy quarks energy 
loss is close in magnitude and thickness dependence to the incoherent linear Bethe-Heithler like limit. This is in 
contrast to the quadratic (EDMS) thickness dependence characteristic for light quarks in the deep LPM regime. 
We note that ref. j34| reached similar conclusions using a different approach. 

At LHC the partial LPM reduction of charm quark energy loss below the incoherent limit is predicted for 
E = 30 GeV jets. However, the predictions for LHC are only quantitatively, but not qualitatively different from 
those at RHIC as long as the densities, and hence plasma parameters at LHC are not much more than about 
four times that at RHIC. 

Finally we show in Fig. 9, the predicted radiation length (see Eq. 114|l ^ as a function of the heavy quark 
mass under estimated RHIC and LHC quark gluon plasma conditions. The main feature to note is that the 
radiation length is comparable or greater than nuclear radii in all cases. The results imply that jet quenching 
is dominantly a volume emission rather than a surface emission phenomenon for all mass jets. 
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For light parton jets, the observed Raa ~ 0.2 accidentally suggest that only jets produced near the surface may 
survive to the detector. However, as shown in ret.) 12). the finite GLV energy loss for jets propagating through 
the entire volume of the plasma quantitatively explains not only the magnitude and pr dependence of Raa{pt) 
in central Au + Au but also its centrality dependence jSH]. Heavy quark tomography will give an important 
complementary test of this physical picture. The computation presented here provides the basis for testing this 
theory on future experiments. The applications of heavy quark tomography to RHIC and LHC experiments is 
the subject of our current work. 
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A Apendix: Assumptions 

To compute the medium induced energy loss, we use the same assumptions as in 16 . For convenience, we recall 
here those assumptions. 

First, we consider a Yukawa potential as in Eqs. ^ and assume that all the Xj are distributed with the 
same density 

p(x) ^ ^p{z) , (28) 

where J dzp(z) — 1. 

Second, we assume that the observed p = [E^, E~ , 0] is high as compared to the potential screening scale, 

i.e. 
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(29) 



We also assume that the distance between the source and scattering centers are large compared to the 
interaction range 



- zq > . (30) 

Finally, we assume that the source current or packet J{p) varies slowly over the range of momentum transfers 
supplied by the potential. 

A major simplification occurs if the relative transverse coordinate (impact parameter) b = — xq varies 
over a large transverse area, A±, relative to the interaction area In this case, the ensemble average over 

the scattering center location reduces to an impact parameter average as follows 

<.->^/^^^... 

The ensemble average over the phase factor then yields 



^^-.(q-q')-b)= (32) 

Also, to calculate the diagrams with one and two scattering centers we need to find the mass correction for 
the full triple gluon vertices A^, and Aij given in However, it is easy to show that in the approximation 
M'^/pI <C 1, these vertices remain the same. Therefore, the full triple gluon vertices including coupling and 
color algebra for producing a final color c from initial color a followed by color potential interactions a„ and 
are then given by 



« -2i.g,i;+fc+e- (k-qi-q2)[[c,a„],a™](T,„(n)T,„(m)) . (33) 

B Diagrams Mi,o,o, ^1,1,0 and Mi,o,i 

In this appendix we present explicit calculation of the diagrams shown in Fig. 8. 



A„ 



A. 
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^ z 1 p ^ |z p ^ ^ P 



S qi'^i I qbai t qi,ai 




k,c 



FIG 8. Three "direct" terms Mi, 0,0, Mi^ifi and Mi, 0,1 contribute to the soft gluon radiation amplitude to 
first order in opacity L/Xoc aei/A±. 

Computation of Mi_o,i 

As a first application, consider the one rescattering amplitude Mi, 0,1- 

Mi,o,i = J ^ iJ{p + k- gi)e^(f+'=-«^)-« Ai(p, k, qi)V{q,)e'^^^^ x 
X iAM{p + k - qi){-i)Amg{k - qi) 
« J{p + fc)e^(^+'=)^° [c, ai]T„, 1 ^ e-*'i-(-^--°) 25, e • (k - qi) x 

X 2E y qi)AM(p + k- q,)A^^ {k-qi) e-^9i-(-i-o) . (34) 

The longitudinal momentum transfer integral, 

/i(p,/c,qi,Zi-Zo)^ J ^v{qi,,qi)AM{p + k-qi)Ara^{k-qi)e-^'''^^''-'"^ , (35) 
can be performed via closing the contour below the real axis since zi > zq. 

In addition to the potential singularities at ±i/ii (/x? = = q? + /x^), the two propagators have two poles 
in the lower qi^ plane, which are approximately located at 

qi = -LUo -Um-ie. , 

(36) 

q2 = -LOo + wi - ze . 

where wq = = ^^''^ ^^'^ '^m = Note that ujm is the term which has the information on 

finite mass corrections. 

The residues give 

gi(u;o +Wm ){zi-zo) 



Res(gi) -v{-uJo - i^m, Qi) 



E+k+{(Ji+U}m) 



gi(wo-tiJi)(zi-zo) 

Res(g2) ~ v{uji - uJo,qi) — — — , (37) 

I^j^ K^[U!i + Urn) 
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while 



" (^2^Ml)i?+fc+(-^m)^ ' ^^^^ 

where we assumed that ^ /ii ^ uji. 

In the well-separated case where iJ,{zi — zq) — 1, this potential residue is exponentially suppressed and 

therefore 



fc,qi,zi - zo) 



X (v{-cJo - Cj^, q^)eH-o+-™)(.i-.o) _ _ q^)g,:(c.o--i)(.i-.o)^ 
t;(0,qi) ^ ' ^ / »(a;o+i),„)(2i-2o) _g»(c^a-^i)(2i-2o)A ^ (39) 



Using the fact that 



2w 
2w 



_ ^ qi)' (40) 



we finally get 



(27r)2 ^ (k - qi )2 + A/2a;2 + ^2 

X (e2b(k^+™^+A^'-')(-i--o) „ g5^(k^-(k-q,)^)(.i-.o)) [c^a^jr^^ ^ (41) 

where bi = xi — xq. 

The Eq. (|41ll represents the massive corrections of Eq. ( A7) in JS] . 

Computation of Mi q o 

Using the same technique as in the previous section we can now compute Afi^o.o- 

^^1.0,0 = / ^JkP + k- gi)e*(P+'=-9i)-« {^g^)e^(^2p - 2gi + fc)" x 



. (2^)^ 

X iAMip -qi+ k)tAM{p ~ qi){2p - qxfV {qi)e'^^''^Ta,aiC 

J(p + fc)e^(P+'=)^°(-*5.aicT,j2i? j ^e-'^^'^U^, (42) 



where 
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Since zi > zq, integral I2 can be performed via closing the contour below the real axis. Only roots {—uq 

2 

: — ie) and {-^ — ie) contribute to this integral (: 
relations given in the Eq. H4n|l this integral becomes: 



ojm — ie) and {-^ — ie) contribute to this integral (root — i/i is suppressed, since ^j,{zi — zq) ^ 1). Using the 



hip,k,ci,,z, - zo) = ^ i,.+i'.+M^,. ^(0^qi)(e^^'^^-^"'-'^"^^^^'"-"^ - 1) (44) 
Finally Mi_o,o becomes 



A/1,0,0 = J(p)e^(^+'=)^° H) / ^v{0,ci^)e-^^^ '^^{2igs) x 

' • ^ (g^(k^W,+M^.^)(.i-.o) _ i)a,cTa,. (45) 



Computation of Mi 1 







Mi,i,o = J -^iJ{p + k~qi)e'^P+''-^'^-»i2p + 2k-q,)° X 

X iAm{p -qi + k)iAM{p + k){igs)€^{2p + fc)"l/(qi)e*«^"ir,,cai 



J(p + fc)e'(f+'=)-«(-zg,r,,mi)(2£; + 2c.) / x 



{p-qi+ fc)2 - AP + ie [p + fc)2 - Af 2 + ie 
Since uj <^ E =^ (2E + 2lj) ^ 2E. 

£c.(2p + fc)" ^ 2pe 
{p + fc)2 - M2 + ie ~ 2pk + k? 

In the large pz limit this becomes 



e^{2p + kr (46) 



(47) 



(49) 



[p + kf-AP+ie ' 4p2 ^ k2 + m2 + Af2a;2 k2 + a;2M2 + m2 

Using these two results, we get 

Mia.o ^ J(p)e'^-°(-»g.)Ta.cai4i^ ^^^ ^ J^' | ^e--^b^/3(;,, fc, q^, ~ zo) 

where 

^3(^'^'^---^") - /t^ (,^,, + fc)2-M2 + .e --'"^^^^-^"^-(-^-^ 

pj(a;o+<im)(2l-2o) 

« (-z) 1.(0, qi) (50) 

2pz 
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Therefore, 



which in the massless hmit, M = nig — 0, leads to Eq. (A8) from |16| . 

C Diagram M2,o,3 

Consider next the gluon two-scattering amplitude JVh^o.s- Fig- 9 shows that for inclusive processes two interesting 
cases arise. 



2,0,3 

Zq Z Zj 



Z2 



^ P 

qi,ai® ® q2,a2 



M' 



2,0,3 

Zq Z Zj 



^ P 

'1.-,™^'! k,c 



qbai< 



q2,a2 



FIG 9. M2,o,3 "direct" contributes to second order in opacity oc {aei/A±)'^, whereas MI^q.s ~ -^^2, 0,3(22 
zi) "contact-limit" contribute to first order in opacity oc (uei/Ax)^. 



In the Feynman diagram approach 



M2.0- 



(27r)4 (27r)4 

X Ai2(p, fc, 91,92) iAM{p + k - 91 - 92)(-«)A„^(/c - 91 - 92)(-i)A„g(fc - 92) 
J(p + fc)e^(P+'^)-'' [[c,a2],ai](T„,(2)r,,(l)) x 



(2^)2 
dqiz dq2z 



(-0 1 ^2*g.e.(k-qi-q2)e-*'i-''^e-''i-''^ x 

(4Sa))i>(qi)w(q2)e-*'?i-(^i~^«)e-*«^-(^2-^«) 



(52) 



(27r) (27r) ((p + fc - 91 - 92)^ + AP + ie){{k - 91 - 92)^ + ml + ie){{k - 92)^ + ml + iej 

where = x— Xq are transverse impact parameters, and we used the soft gluon and rescattering kinematical 
simplifications Eqs. ll^ . e.g. J{p + fc — 91 — 92) ~ J{p + k) fa J{p)- For the qiz integral, it is convenient to 
rewrite the phase as 

g-i9lj(2l-2o)g-j92z(z2-zo) _ g-j(9lz-|-922)(zi-2o)g-i92j(z2-zi) 
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The first longitudinal integral is closely related to Eq. H35|l 



Since 21 — zo ^ 1//^: we again close the contour in the lower half qiz plane and neglect the pole at — i/ii. 
The remaining qiz poles are shifted by —q2z and qi ^ qi + q2 relative to Eq. (|^ : 

gi = -q2z - Wo - - ie , 

(54) 

q2 = -q2z - t^o + W(i2) - «e , 
where W(i2) = ^ ^^2uj^^ ■ '^^^ residues at qi, q2 then give 

I {v{-q2z ~ Wo - w™, qi)e»('^o+i'™)(^i-3o) „ „(„^2. - wo + W(i2), qi)e'('^«-'^(i2))(^i-^«)) 

^2 ~ 77+; — \ ' (.55) 

where we have neglected C'(exp(— /iA)) contributions. This differs from Eq. H39|l mainly in that the potential 
is evaluated near — (?2z, which still remains to be integrated over, and wi — > lo(^i2)- 
Next we need the following critical q2z integral 



dq2z v{-q2z + Suj,qi)v{q2z,(i2)e-"^^''^^^'^'^ 
27r {{k - q2)'^ + m? + ie) 



/3(fc,qi,q2,22"^l) = / — x2 , ^2 , ■ (56) 



In the general case (including the special contact case with Z2 = zi) both q2z = —ilJ-2- ~*Mi singularities in 
the Yukawa potential contribute together with the pole at q2z — uj2 — luq — ie, resulting in 



/3(fc,qi,q2,Z2 - zi) w — u(0,qi)w(0,q2) e 



(47rQ!s)2 /g-A'2(z2-Zl) g-/Jl(22-Zl)g-il5w(z2-Zl) 



(57) 



Fortunately, we are interested in only two extreme limits: 

• The limit of well-separated scattering centers Z2 — zi 3> 1//^ ; 

• The special "contact" Z2 = zi limit to compute unitary contributions. 
For Z2 — zi = \ 1 / ^ the Eq. (|57|l reduces to 



/3(fc,qi,q2,Z2-^i»l^) « --^z;(0,qi)z;(0,q2)e-'("^-'^°)(-^--i) . (58) 



For the special contact contribution Z2 — zi = it reduces to 



/3(/c,qi,q2,0) « _L„(0,qi)«(0,q2) . (59) 



i.e., exactly ^ of the strength in Eq. (|58|l . 
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The contact limit of this amplitude is therefore 



^.(0,q,)e ^H)y 
Note that in the massless limit it reduces to Eq. (BIO) from 

D Diagrams M2 0,0 and M2 2,0 

In those graphs it is the jet rather than the gluon that suffers two sequential scatterings as seen from Fig 
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2,0,0 ^< 



k,c 



qi.ai 



Z2 



q2,a2 



M 



2,0,0 



k,c 
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q2,a2 
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Zl 



k,c 



f 2 



qi.ai 



I) ^ 

^ q2,a2 




FIG 10. M2,o,o and M2,2,o graphs in the well-separated case together with their 22 ~ z\ limits M^^fi, 

M2.o,o = f fy^, iJip + k-q,- g2)e'(P+^"^^-^^)"°l^(gi)e'^^^^ V(g2)e^^^"^ x 

J [ZttP [ZttP 

X {-i2Efigs{2p + k) iAMip + k - qi - q2)iAM{p - qi - q2)iA.M[p- q2)a2aic{Ta^Ta^) 
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dqi^ dq2z u(qi)w(q2)e-^'!'i.(^i-^o)g-i«2.(z2-^o) 



(2^) (2^) ((p + A: - 91 - 92)2 + + _ _ q^)2 + m2 + ze)((p - 92)2 + KP + zej 

In this case we define 

/2(p, fc, -1 - -0) = y ^ ^ ^ _ _ ^^^2 + + ,,)((^ _ _ + + • (62) 

Since zi — zq ^ 1//^, we neglect the pole at — i/ii. The remaining qiz poles are 

(63) 

(72 = -922 - «e , 

where we discarded (p + k — qi — q2)^/E'^ relative to luq. The qi, 92 residues then give 

w(-g22,qi) |^g»(a;o+i:.,„)(2i-2o) _ 1^ _ (64) 



l2~-i 



Not that Wo has been neglected in the potential relative to /ii. In the contact limit, the second integral, I3, 
is then equal to 



hiP,<ll,<i2,Z2- zi) = J 



dq2z ^(-922, qi) v{q2z,ci2) e-"?^-(^^-^i) 
(27r) ((p-g2)2 + Af2 +ie) 

I I 1 if uA = /i(z2 — 2i) ^ 00 



2 

With the help of Eqs. (|64I65() in the case of contact limit we obtain 



M2,o,o = ^J(p)e'(^+'=)- / ^^.(0,qi)«(0,q2)e-'('i^+'i^)-^^ x 



2 '-^^ y (27r)2 (27r 



In the same way we will obtain for M2,2,o 



M2^,2,o = ^^b)e^^^+'=^^°/^^-(0,qi)z;(0,q2)e-(''^+'>^)-^^ X 



" ('^g-) k2+i'2+^M2x2 e-(''^+™^-^^^)(--) ca2a, (r.,T.J . (67) 
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E Diagrams M2 0,1 and M2 0,2 



Bellow we compute the case when one of the hits is on the parent parton and the other hit is on the radiated 
gluon. Explicit calculation is shown on the example of M|g i — M2^o. 1(2:2 — zi) in Fig. 11. 
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qi'^i ^jI IMI J^ q2,a2 



FIG 11. Af| o 1 and M| o,2 topologically indistinct contact diagrams. There are no additional factors of | 
arising from the integration in taking the Z2 — z\ limit in Af2,o,i and iV'/2,o,2. 



X (-i£:+)Ai iAM(]3 + fc-9i - 92) (-i)^mg(fc- 9i)«^Af(p- 92) 



2 X 



X 2i5,(e- (k-qi)) e^"«^° a2[c, ai] (T,,r„ J (2£;) 

I dqi^dm^ «(qi)i>(q2)e-»g^-(^i-^°)e-'^^'(^^-^°) 

/ (2^) (27r) ((p + fc-gi -q2)2+M2 + ze)((fc-gi)2+m2+ie)((p-92)2 + M2+ze)-^ ^ 



We perform the q\z integral first 



/2 fc,qi,q2,zi - zo) = / ^^77 — X2 I A^2 I ■ — 2 I ■ ^ ■ ^9 

The pole at — i/ii is again exponentially suppressed. The poles of interest in the lower half plane are 

q_\z = -q2z -ujQ-UJm-ie and qi^ = -ujq + uji - ie. 
Taking the residues leaves us with 

I2{p,k,qi,q2,zi - zo) = 



E+k+{q2z +uji +u;„i) 

X (t;(-g2.-c^o-<:^m,qi)e'("°+^")("^-"")-t;K-c^o,qi)e^^"«"''^^""^^'^^"'«^) ■ (70) 
It is important to notice that there is no pole at q2z — —^1 in Eq. (|70|) . The remaining integral over q2z is 
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^3(p, ^,qi,q2,zi - 20,22 - 2i) = J 



Stt q2z +UJ1+ oj„i ^ {p - (72)^ + + ie 



Q-i{q2ziz2-Zo) + {uJl-UJo)(zi-Zo)) 

X v{-q2z - uJo-Um,q.i)v{q2z,q.2)) 7 „ , .,0 , . v{uji -uJo,qi)v{q2z,q2))- (71) 

The poles in the lower half plane are q2z = — ie, q2z — and q2z — —ojq — Qm — il^i- In the well 

separated case, contributions from second and third residue are exponentially suppressed oc exp[— /i2(zi — zo)], 
and therefore can be neglected. In the contact limit contributions from second and third residue cancel exactly 
(for more details see Eq. (d5) ^^l- Therefore, we get 

Ml,^, = J{p)e'^^+^^^'>{-i) I ^e-'>^-^^«(0,qi)(-*) / ^e-'^-''^^;(0, q2) x 
x2^5.71 J. ..ri , , {e'' ^ ?2.)(..-.o) _ a2[c,a,]{Ta,TaJ .(72) 



(k - qi)2 + M^x^ + 
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In the massless limit this equation reduces to Eq. (D6) from |16| . 

Notice that, unlike to the previous examples, there is no factor of ^ in M2 q i- 

We can get M|q 2 from q by replacing every 2 with 1 and vice verse. 

However (since scattering centers are identical), we need to symmetrize this two diagrams, which effectively 
leads to multiplying every diagram with ^ . 

F Zero measure contact limit of M2 i o and M2 1 1 

In calculating the different contributions coming from two interactions with the same potential centered around 
Xi we have to take into account the two graphs given in Fig. 12, where one of the hits occurs before the gluon 
emission vertex and the other one after. 




M2,u 



qbai 




^ p 

J k,c 
q2,a2 



FIG 12. Diagrams M2^i q and M| j i of 0(0) according to the time-ordered perturbation theory. 

In the framework of time-ordered perturbation theory of li-jyilHSj the graphs are identically zero in the contact 
limit because J^^ dt--- = 0. Here we present a more detailed study of the validity of this argument in the case 
of massive quarks and gluons. 
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(27r)4 (27r)4 

X iAM(p + - - g2)l^(gi)e^«^^^iAM(p + fc - q2){2p - 2q2 + k^e^iAMip - 92)i5«'5^(?2)e'«=^= x 
X {-i{2p+2k -qi- 2q2f){-i{2p - g2)°)a2caiT„,T„, 
J{p + fc)e^(f+'=)-''a2cai(T„,T„J(-i5«)(2^)2 x 



/ ^^^(^i)«"'^'^'''"''°^^(^2)e-'^^^''=-''°H2f'-2(Z2 + fc)'^e^ x 



X Am(p + fc - gi - q2)AM{p + k- g2)AM(p - 92) 



(27r)2 (27r) 



— (2p - 2(Z2 + kYe^ _ ^^^2 _ ^2 + + _ ^2)2 _ M2 + le) J 27r ((p + fc - g)2 - + ie) 

= J(p + fc)e'(?'+'=)-°a2cai(T„,T„J(-i5,)(2E)2y'^^e-''>^''^e-''>=''=^ x 
dg2z v(q2)e-'9^'(^^-^i) 



/ 



27r ((p - 92)2 _ M2 + ie){{p + k- 92)2 - M2 + ie) 



x-.( 1- ,2.q0e^^^^--^ (73) 

2 

In the contact limit, there are four roots which contribute to the integral over (— ^^2? — i^i, tt^ and 
(-Wo - <^m)), leading to 

J2, 



M|,i,o = J(p + fc)e^(^'+'=)-''a2cai(T„,T„J(i5.)y'^^e-'('>^+''^)b^ 



X t;(0, q2)t;(0, qi ) -^e' 



k2+m2+M2x2 2i 2oj MiM2(a<i + At2) 



~ ^J(p+fc)e^(f+fc)-oa2cai(r„,T„J(i5,) x 
M 

7(277)2(277)2'^ ^ '^'^^ '^'^k2+m2 + M2a;2^ ^'^^ 

Therefore, the contributfon of M2_j q is suppressed by an O(^) factor relative to the other graphs. 



X iAM(p + fc-gi -(72)(-i(2p + 2fc-(7i -2(72)")y(<Zi)e'«^"^zAM(p + fc-g2)(-iA„^(fc-g2)) x 
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J{p + fc)e*(P+'=)-" {-2ig,){2Ef[c, a^^Ta^Ta, x 
X Aa/(p + k- qi- q2)AM{p + k~ q2)A.,n^{k - (72) 

J{p + k)e^(^+^'>^»{-2igs){2Ef[c,a2WTa,Ta, j ^^e-^^-'^-e-^^^'^^{e ■ {\^- c^)) x 
dq2z t;(q2)e-'«^^(^^-^i) Z" dq, v{q, - (?2., qOe'^'^^^^-^") 



27r {{k - q2y - ml + ie){{p + k - q2)^ ~ + ie) J 2tt {{p + k - q)^ - A'P + ie) 
= J(p + fc)e^(P+'=)-°(-2*5s)(2i?)2[c,a2]aiT„,T,, J ^^e-<iibie-<i^b.(^ . _ q^)) x 
dg2. t;(q2)e-*«^-(^^-^i) 



2tt {{k - q2Y ~ m?g + ie){{p + k- 92)^ - M'^ + ie) 

2^"^^ fc fe, qi)e' (^^-^"^ (75) 

2 

In the contact limit, there are also four roots which contribute to the integral over 522 {-~ip2, ~*Mi7 ^f" ^-i^d 
(ui2 — Wo)). After computing the residues, we get 

^^2,1,1 = Jb + A;)e'(P+'='""(2z.g,)[c,a2]aiT,,T,, x 

/■ d^qi d^q2 (^,+^,)b, (£-(k-qi) 
(27r)2(27r)2' "^^^k - qi)2 + + M^x^ 



2i 2a; MiM2(mi + M2) 



X e 



i^J(p + fc)e*(P+'=)"°(2ig,)[c,a2]aiT,,Ta 

/i 



(2'kY (2itY k^ + m^ + Ai^x^ 

We see that the contribution of M|]^ p is also suppressed by an O(^) factor relative to the other graphs. 
Therefore, we will neglect the contribution of M| ^ ^ and 1 in the energy loss calculation. 

G Computation of the first order radiative energy loss 

In this section we want to compute the first order in opacity radiative energy loss. According to Eq. (|10|l we 
have 

SNf^VN, = i(T,(|M,P> + l.ReT.(M„.Ay,^^, ,77) 



28 



Ml is sum of all diagrams with one scattering center and M2 is sum of all diagrams with two scattering 
centers. 

Using the results for Afi_o,o, Mi^i^o and Mi^o,i from Appendix A we get 



Ml = Mi,o,o + Mi,i,o + Mi,o,i = J(p)e^(f+'=)^°(-*)(2»5.)T^aiy ^v{0,qi)e-"^'''' 

X {( (g-(k-qi)) <^ 'ii ^ i(u„+u,rn)(zi-z„)i 1 _ 

^^(k-qi)2+m2+M2x2 k2 + m2+M2x2^ ^ ' 

_ (e- (k-qi)) i(a,o-c.i)(^i-^o)rc a,l llif a^cl f78l 

(k-qi)2 + m2+M2x2'^ ^'''"'^^ k2 + m2+M2:c2«iC^ 

where 

Wl - ' ^™ - 



Then, 



^(iMip) = iV|J(p)p(45.^)J-|^Kqi)P 



C2(T) 



rg,/ (e-(k-qi)) .2 , g / (e-(k-qi)) 2 

^ ^(k-qi)2 +m2 +M2x2 k2 + m2 + M2a;2 ^ ^ {k - q^)^ + + M^x"^ ' 

(e-(k— qi)) e-k (e-(k — qi)) n tt - \( 

(k - qi)2 + m2 + M2a;2 ~ k2 + m2 + M2a;2 ' (k - qi)2 + ^2 + M2a;2 + '^"^^^^ ~ 

~ (k - qi)2 + m2^+ M2x2 ~ k2 + ml ^ M2x2 ^ k2 + m2 ^ M2a;2 ^ + ' ''^^ + 

+ " k2 + ml + M2.2 (k - q02 +;n2 + M2.2 ' - "^^(^^ - + 

+ ^k2+m2!M2a.2)'^™'^'> 

Here we used Tr(TajTo2) = ^^^^^^^^^ , and defined a = Tr{(?a? — cam). Factor A'' comes from sum over 

A^ scattering centers. 

To compute M2 we will first add all the diagrams with two scattering centers calculated in Appendix B-E, 
and than take their average. 



M2 = ^J(rte'(''+*^)"'(-2*5«)Ta,T,,-^ 
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X ^ k2 + J2 + M2x2 {e'^"°+"-^<^^"^°^([[c, a2], ai] + [a2ai,c]) - [[c,a2],ai] - asaic} - 
" (k - qi)2 + m2 + M2x2 " e'^-o-..)(..-.o))^^^^,^ + ^^[,^ ^^])} (80) 

Using this we can now find ^ReTr(MoM2). 
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^-ReTr{M*M,) = 7V|J(p)p(4a-^ /^Kqi)P^'^^^ 



2^. e-k (e-(k-qi)) 



k2 + m2 + APx^ (k - qi)2 + m2 + APx^ 

X {COS((WO +il'm)(2;i - ^o)) - COS((WO - CJi)(zi - Zo))}} (81) 

Therefore, ^ (|A/i|2) + ^ReTr (M^Ah) is equal to 

^(iMiH + AReTr(M*M2) = Z?«|J(p)P(4a^C^^|^Kqi)P 

r , (g-(k-qi)) e-k (e-(k-qi)) 



(k - qi)2 + m2 + APx'^ ^k^ + + Af2^2 (j^ _ q^)2 _|_ ^^2 ^ ^22,2 ' 

(k - qi)2 + ml + APx^ 
X (l-cos(^ 2 (zi-zo)))}. (82) 



Here we nave used that wi + uj,n — 2^7^ and a — ■^Cj^Un 

Using Eqs. I^|H1) it is now easy to extract dE^^^ = tod^Ng 



rr,r^a L ^ r d^qi ^2 M^k ^ 



dx TT X J TT (q2 + ^2)2 



TT 



r o (e-(k-qi)) . e-k (e-(k-qi)) 



(k - qi)2 + to2 + APx^ + m2 + M2a;2 (k - qi)2 + m2 + M2a;2 

/• (k-qi)2 + m2+Af2a;2 e"^^^ 
X / dz,{l - cosC- 2 (^^ _ ^ , (83) 

where we have used Eqs. (1, 6) from |16| to write the resuh in terms of opacity L/A. For simpUcity we have 
also assumed exponential distribution exp(— Az/L)/L between scattering centers. After the z\ integration we 
get 



d4nd _ CnasL^ f dk^ f d^m ^^2 



dx TT A 7 k2 + to2 + M2a;2 J tt (qf + fj.'^y 

2 k ■ qi (k - qi)2 + {ml + M2a;2)qi ■ (qi - k) 

(^)' + ((k-qi)2 + M2x2+m2)2 ^ ' 



which in the massless limit reduces to Eq (125) from [16|. 
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